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For discrete-time systems, governed by Kraus maps, 
the work of D. Petz has characterized the set of 
universal contraction metrics. In the present pa- 
per, we use this characterization to derive a set of 
quadratic Lyapunov functions for continuous-time 
systems, governed by Lindblad differential equations, 
that have a steady-state with full rank. An extremity 
of this set is given by the Bures metric, for which the 
quadratic Lyapunov function is obtained by inverting 
a Sylvester equation. We illustrate the method by 
providing a strict Lyapunov function for a Lindblad 
equation designed to stabilize a quantum electrody- 
namic "cat" state by reservoir engineering. In fact 
we prove that any Lindblad equation on the Hilbert 
space of the (truncated) harmonic oscillator, which 
has a full-rank equilibrium and which has, among its 
dccoherence channels, a channel corresponding to the 
photon loss operator, globally converges to that equi- 
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1 Introduction 

The study of asymptotic convergence properties of 
open quantum systems has been reinvestigated in 
the last decades thanks to the development of quan- 
tum information [12] and reservoir engineering [15] 
to stabilize and protect fragile quantum states (see 
for example [7] and the references herein). The 
most standard models assume memoryless environ- 
ments and are described by Lindblad differential 
equations (continuous-time) or Kraus maps (discrete- 
time). During the 70's, a line of work around 
[3, 20, 5] has developed a sufficient algebraic con- 
vergence criterion for these models, known as the 
Davies-Sphon-Frigcrio (DSF) criterion. It ensures 
existence, uniqueness and global attractiveness of a 
full rank steady-state density operator. Meanwhile, 
the particular structure of the quantum master equa- 
tions also allows to investigate convergence with a 
geometric approach. Indeed, one can define "univer- 
sal" contraction metrics which are independent of the 
precise dynamics [13]. As noticed in [11], contrac- 
tion analysis opens the way to developing systematic 
convergence-characterizing tools by examining where 
the contraction is strict . 

In the present paper, we propose such a tool in 
the form of a systematic Lyapunov function design 
for quantum systems described by Lindblad differ- 
ential equations; and we illustrate it by applying it 
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to a practical reservoir engineering situation, moti- 
vated by our recent reservoir engineering proposal in 
[18, 17]. We limit ourselves to finite dimension for 
simplicity, although in principle similar results should 
be possible in infinite dimension. We thereby start 
from a characterization of the metrics that are (non- 
strictly) contractive for all Kraus maps [13]. This 
characterization then leads to a systematic construc- 
tion of (non-strict) Lyapunov functions for all Lind- 
blad equations with a full-rank equilibrium; a main 
advantage of these Lyapunov functions is that their 
time-derivative has an explicit expression which lends 
itself well to an analysis of critical points. The result- 
ing convergence criterion is different from the DSF 
criterion since it does not assume that the vector 
space spanned by the Lindblad operators is closed 
by Hermitian conjugation (modulo identity). This 
allows for instance to rapidly show the following re- 
sult not covered by the DSF criterion, in the context 
of our engineered reservoir application: any Lindblad 
equation on the Hilbert space of the (truncated) har- 
monic oscillator, which has a full-rank equilibrium 
and a decoherence channel proportional to the pho- 
ton loss operator, globally converges to that equilib- 
rium. 



The paper is organized as follows. Section 2 
recalls the models describing the evolution of open 
quantum systems, both with discrete-time and 
continuous-time dissipative channels. Section 3 
reviews known contractive metrics for discrete-time 
quantum systems, described by Kraus maps. Sec- 
tion 4 contains our main results about systematic 
Lyapunov functions for continuous-time quantum 
systems, described by Lindblad-Kossakowski dif- 
ferential equations. Section 5 applies the result to 
obtain a convergence proof for a reservoir engineered 
to stabilize quantum cat states of an electromagnetic 
field mode. 



2 Dynamics of open quantum 
systems 

All the material presented in this section can be found 
with much more details in [6] with developments cen- 
tered around the experiment described in Section 5 
and in [1, 21, 2] for more formal theoretical exposures. 

We consider a quantum system on a finite- 
dimensional Hilbert space H, whose state is repre- 
sented by a density operator p, i.e. a self-adjoint, 
unit-trace positive semi-definite linear operator on % 
which represents quantum probabilities. We denote 
the set of all density operators on ri by V-^. If the 
quantum system is closed, it evolves according to the 
Hamiltonian dynamics (Schrodinger dynamics with 
h=l) 

where i = s/^l, [A, B] = AB - BA and the Hamilto- 
nian H is a self- adjoint operator on %. Such evolution 
leads to a unitary propagator, i.e. the evolution is an 
isometry for each t and two states never converge to- 
wards each other. 

This behavior can be changed by considering open 
quantum systems. In particular a measurement op- 
eration on a quantum system perturbs it in a non- 
unitary (in fact stochastic) way. The present paper 
in contrast considers open systems where information 
is lost to the environment without being measured. 
The latter can be conceived as the reduction to our 
target system, on Ji, of a Hamiltonian dynamics on a 
larger Hilbert space T-L =H<£) H which is the tensor 
product of the target system and an external "bath" . 
Open quantum systems can both be helpful, for sta- 
bilization, and detrimental if the bath induces dissi- 
pation "in the wrong direction" w.r.t. our goal, also 
known as decoherence. Typically, a quantum system 
interacts with both undesired decoherence sources, 
and stabilizing baths designed to counter them. 

One of the mainstream open quantum models as- 
sumes that the environment loses any information 
about the system between consecutive interactions. 
In discrete-time, this Markov assumption leads to 
a trace-preserving completely positive map or Kraus 
map $ for the density operator of the target system 
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[9]: 

p -> $(p) = ^M,pM] (1) 

ft 

where ^ fc A/^/ fc = I . (2) 

The Mfc are arbitrary linear operators on H (but sat- 
isfying condition (2)), t denotes operator adjoint 
and I is the identity operator. In continuous-time, 
the Markov assumption leads to a class of models 
known as Lindblad(-Kossakowski) differential equa- 
tions [10, 8]: 

TtP = C iP) = ~*W,p] ( 3 ) 

-fE(4^+p4i fe -2L feP 4). 

k 

The I/fe are arbitrary linear operators on H. In the 
following sections, we review contraction-like tools to 
assess the convergence of (1) and develop how they 
lead to Lyapunov convergence arguments for (3). The 
stationary points of such dissipative evolutions are 
also known as pointer states in the physics literature 
[22]. 

A particular case of Markovian open quantum sys- 
tems is obtained through reservoir engineering. Wc 
start with a discrete reservoir. In this case, at each 
step, a 'new' auxiliary system (thus carrying no mem- 
ory about past system states) is prepared in a pre- 
determined initial state p B {0) and brought into con- 
trolled Hamiltonian interaction with the target sys- 
tem for a time T; the coupled systems thus evolve 
as 

p®p B (0) -> U T {p®p B {$))U ] T 

where Ut = U(T) is a unitary operator onH® H , 
solution at time T of the controlled Schrodinger equa- 
tion 

-^U(t) = —i H controlled (*) U 

with initial state U(0) = I. Target and auxiliary sys- 
tem are in general entangled. After interaction the 
auxiliary system is discarded; the state of the tar- 
get system is then described by taking the partial 
trace of U T (p <£> p B {0)) U\. over the auxiliary system 
Hilbert space, and the resulting effect is described by 
a Kraus map of type (1), where the Mk depend on 



p B (0) and the function i?controiied(i) over [0,T]. Now 
a new auxiliary system and interaction can be initi- 
ated, to iterate the (potentially time-varying) Kraus 
map. With proper tuning, a stabilizing effect can be 
obtained (see Section 5), without using measurement- 
based feedback. 

In the limit of a very large number of short 
interactions per time unit (T = di < 1 and U(T) 
close to I), a continuous-time model of type (3) can 
be obtained. The latter is also often viewed as the 
continuous interaction of the target system with one 
large infinite-dimensional system (e.g. thermal field 
environment). 

3 Contractive metrics for 
Kraus maps 

The content of this section is also presented with 
much more details in [7] and relies essentially on the 
key contributions [14, 13]. 

Definition III. 1: A distance function d : pi,p% — > 
d(pi, P2) G R defines a contractive metric for a Kraus 
map i> if and only if d($(pi), $(p 2 )) < d(p 1 ,p 2 ) for 
all pi,P2- 

Motivated by popular quantum results, we first 
consider 'distance functions' in a wide sense: they 
may be asymmetric with respect to their two argu- 
ments, sometimes strictly defined only for full rank 
pi and/or p%. The latter point may be annoying 
when the goal of an evolution is to reach a pure state, 
i.e. a rank one density operator; in presence of dis- 
turbances, a small residual proportional to I usually 
solves the situation. 

Thanks to the particular structure of Markovian 
open quantum systems, there exist several metrics 
which are contractive for any Kraus map; see e.g. [7] 
for a summary: 

• The trace distance dtr{pi,P2) = 
lTr(V(pi-p 2 ) 2 ). 

• The Bures distance dB(pi,P2) = 
y/l - F(pi,p 2 ), where F(p 1 ,p 2 ) = 
Tr ( y/y/ptpi \fpi) is the quantum fidelity. 
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• The Chcrnoff distance dc{p\,P2) = 
i/l - Q{pi,p 2 ), where Q{pi,p 2 ) = 
min < s <i Tr(pf p\~ s ) . 

• The quantum relative entropy ds{pi,P2) = 
v'TrO^logpi - logp 2 ))- 

• The quantum x 2 -divergence d x -2.(pi,p 2 ) = 
1^Tr( (pi - p 2 )P2 1/2 (pi - P2)P2 1/2 )■ 

• Hilbcrt's projective cone metric: dn(pi,P2) = 

1 I \ I -V2 -l/2\ / \ / -1/2 -1/2 U 

log(A max (p 2 pip 2 )/A min (p 2 Pl p 2 )) 
if pi and p 2 have the same support, and 
otherwise dn(pi, p 2 ) =00. 

All these metrics have their advantages and draw- 
backs, and none of them is always strictly contract- 
ing. The Hilbert metric for instance is symmetric 
but non-Riemannian, and also contractive in the dual 
(a.k.a. Heisenberg picture): defining 

&{X) = Y,MlXM k (4) 

k 

for each self- adjoint operator X on H, such that 
Tr($ * (X) p) = Tr(A $(/?)), we have that $*(I) = I 
and that dn is contractive for $*: this shows 
that any $* tends to bring the eigenvalues of any 
operator closer to each other. Furthermore, the 
contraction ratio for 4> is given by tanh(A($)/4) 
where A($) = max{<i_f/(pi, p 2 ) ■ Pi,p 2 > 0}. Sec 
[19, 16] for recent illustrations. 

When considering deviations from a target solu- 
tion, it can be attractive to compare p and p + Sp 
with a Riemannian metric on the set of strictly pos- 
itive Hermitian operators of trace one. This set is a 
sub-manifold of the vector space of Hermitian oper- 
ators. For each p £ with full rank, Sp belonging 
to the tangent space at p, is any trace-less Hermitian 
operator. By linearity of Definition III. 1 readily 
yields: 

Definition III. 2: A Riemannian distance metric \\ ■ 
\\. : p,Sp £E T>fi x TVfi — > \\5p\\ p G R is contractive 
for $ if and only if (with slight abuse of notation) 
||$(<5p)||$( p ) < ll^pllp for all traceless Hermitian Sp 
and all p > 0. 



The following proposition results from the remark- 
able result by Petz [13, 14]. 

Proposition III. 3: Consider the Riemannian met- 
rics that are contractive for all (finite- dimensional) 
Kraus maps. The set of these metrics can be param- 
eterized via the formula 

\\8pf = f Tr(6p{8uj s + 8u>l)/2)m(s)ds (5) 
Jo 

where 8uj s is solution of the Sylvester equation 

s Su! s p + pSuj s = Sp (6) 

and m s ds is a positive finite measure. 

Proof: In [13] the metrics are parameterized by 
standard operator monotone (decreasing) functions, 
that is functions / on the positive cone satisfying 
f(A) > f(B) if A < B (A, B Hermitian positive 
operators), with xf{x) = f{l/x) and /(l) = 1. 
Such functions can be parameterized by positive fi- 
nite measures m(s) for s € [0, 1] such that 

Introducing right and left multiplication superopera- 
tors L p a = pa and R p o~ = o~p, the metrics in [13] are 
written: 

\\5p\\ 2 p = TtiSpfiLpR-^R-Hp). 

Plugging (7) into this yields the announced form. □ 
The two extreme metrics in this set are as follows. 

• A measure m(s) concentrated on 1, gives the Bu- 
res metric Tr(SpSuj) = 2Tr((5w 2 p) with hermi- 
tian Sui satisfying pSui + Scop = Sp . Its geodesic 
distance is the Bures distance ds- 

• A measure m(s) concentrated on 0, gives the 
symmetric metric Tr( 8pp~ l 8p ). 
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4 Contract ion- based Lyapunov 
functions for Lindblad dy- 
namics 



Our goal now is to make use of these Kraus map con- 
tractive metrics to study convergence of Lindblad- 
Kossakowski differential equations. In particular, we 
want to improve the DFS criterion [3, 20, 5] which 
gives a sufficient algebraic condition characterizing 
the omega limit set. We therefore look for a system- 
atic and strong Lyapunov function design method. 
For standard characterizations of the distance to a 
targe t state ptarge t, like e.g. fidelity F(p,p taxget ) = 
Tr ( y/ y/PPtaxgct \ff>) , the time-derivative can be hard 
to compute in an operator setting. It is therefore 
meaningful to look for Lyapunov function construc- 
tions with handy computational properties. Our 
main results build up as follows. 

Theorem IV. 1: Consider a full-rank density oper- 
ator p subject to any Lindblad equation of the form 
(3). Consider any of the Kraus map contractive met- 
rics parameterized as in Proposition III. 3. Then such 
metric is also contractive under the Lindblad equation 
as, for p and bp following solutions of the Lindblad 
equation, we have for each fixed s: 

iTr(6p6u s ) = iTv(6pSJ s ) 

-sJ2 Tr([5uj,L k }p[6co,L k ?) 

k 

-Y,Ml5^,L k ]p[S^,L k ?) <0 (8) 

k 



with 8u) s (p) solution of (6) at each time. 

Proof: It suffices to establish (8), since then the 
contraction of the metric follows from integrating the 
negative function defined there over m(s)ds. 

Consider the function f s (p,Sp) = Tr(8pSuJs)- We 



have, using (6) and its adjoint, 

£j s (p,8p)=Tr(^8u> s + 8p^) 

= Tr( + to\ (P^ + s^p) ) 

= Tr( q*5u a + 8u\ - - %6u a ) 

= Tr( *§£ (Scj s + Scot)) - % (aSulSu. + SujM) ) ■ 

(9) 

Now the Lindblad equation (3) can be 
plugged in for p and 5 p. The term 

in [H, p] yields a term proportional to 
Tt(B — B^) with 

B = 5pH6u s + 8pH6oj\ - {s6w\8u s + 8uj s 8w\)pH . 

Replacing 8p by 8ui\p + sp8uj\ in the first term and 
by s8u) s p + p8ui s in the second term, one easily sees 
that Tr(B) = Tr(B t ). So the Hamiltonian term has 
no effect on j- t f s {p,8p). 

A similar procedure can be followed for the terms 
in Lfc, for each k. Namely, once (3) is plugged into 
the last line of (9), replace Sp respectively by its ex- 
pression as a function of 8ui a , p in the terms contain- 
ing 8w\ and by the adjoint expression (as a function 
of 8oj\ , p) in the terms containing 8lu s . Then several 
terms simplify, grouping the remaining ones precisely 
yields the right hand side of (8). Since the content 
of the trace operator there is self-adjoint, the same 
result must hold for /J = Tr(SpSujl). □ 

Corollary IV. 2: In particular for the Bures metric, 
we have 

iTr(SpSoj) = -2^2lr([5u,L k ]p[Su,L k ]^ <0 

k 

with pSui + Sup = Sp. 

This result can be used to systematically build a 
Lyapunov function for any Lindblad differential equa- 
tion, for which an equilibrium p^ (and nothing more) 
is known. Namely, we take poo as basis for the tan- 
gent vector (p — poo) in the above results. 

Corollary IV. 3: If the equation (3) has a full-rank 
equilibrium p^ , then whatever the form of H and of 
the L k , the function 
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Figure 1: Simplified scheme of the cavity quan- 
tum electrodynamics experiment at Ecolc Normale 
Superieure, Paris. 

V Bmcs {p) = Tr( Poc G 2 p ) (10) 
with poo G p + Gppoo = p - poo (11) 

is a (non-strict) Lyapunov function for (3) and: 

iv Bmcs (p) = - E Tr ( t G p> L ^ p°° t G P' L *P ) ■ 

k 

(12) 

In order to conclude about global convergence of 
the state to Poo , it then remains to examine how 
the commutator of the set of all Lk relates to poo , if 
necessary with a LaSalle-type argument. 



5 Application: a reservoir for 
quantum cats 

5.1 Physical description 

We apply the above framework to a cavity quantum 
electrodynamics experiment, whose goal is to manip- 
ulate the (quantum) state of an electromagnetic field 
mode through its interaction with atoms. A scheme 
of the setup is shown on Figure 1; see [6] for a thor- 
ough explanation of its working. 

In summary, a microwave field at frequency u> is 
trapped in the cavity C made of superconducting mir- 
rors. Atoms are sent one after the other through C, 
with the goal that the field interacts with an atomic 
transition at frequency ujq w lu. To this end, each 
atom is initially prepared by B in one of the two 
states concerned by the relevant transition, and a 
classical microwave pulse in R\ allows to put it in 



any superposition of those two states, parameterized 
by u. In addition to choosing with u the initial state 
of each atom, we can control the experiment by im- 
posing during each interaction a tailored evolution 
of 5(t) = uio(t) — w, the detuning between field and 
atomic transition frequency, with a good time reso- 
lution (via a Stark effect on the atom). Varying S 
around allows to combine resonant and dispersive 
physical effects, on the way to producing truly quan- 
tum phenomena (see below). 

The elements i?2 and D can be used to detect the 
atomic state e.g. for feedback purposes [4]. In the 
present case they play no role, since we use the atomic 
interactions as an engineered reservoir to control the 
field. This means that we consider the Kraus map as- 
sociated to the field evolution when tracing over the 
atomic states. In addition to this (supposedly) stabi- 
lizing action, we want to consider the effect of deco- 
herence due to spontaneous photon loss of the field, 
characterizing how the engineered reservoir allows to 
counter it. We therefore begin by describing the sit- 
uation without decoherence in discrete-time (Kraus 
map), then derive its associated Lindblad equation, 
and add the traditional decoherence channel to this 
Lindblad equation. From there, we can compute an 
equilibrium and analyze convergence using the result 
of Section 5. 



5.2 Goal; Kraus and Lindblad models 

An electromagnetic field mode at frequency u is mod- 
eled as a quantum harmonic oscillator. Its Hilbert 
space is spanned by the orthonormal basis of Fock 
states \n) with n = 0,1,2,..., which are the eigen- 
states of its closed-system Hamiltonian: 



IT 



hu(n+ 1/2) \n)(n\ =I/2 + SwN. 



n=0 



Here we have defined the photon-number operator 
:0 n\n)(n\ and the identity operator I — 
J2n=o l n )( n l- We will denote the generic field state 
by p or \tp). 

The closest in this setting to a classical field state 
of complex amplitude a = Ae 1 ^, is the so-called co- 



G 



t) = eH"l 2 /2 £ « | n) . 



herent state 



To highlight specifically quantum features, in [18, 17] 
our goal is instead to reach quantum superpositions 
of such coherent states, of the type 



N 

£ 

*;=i 



S k2in/N\ 



(13) 



with fixed and known /3 k g C, |/3 fc | 2 = 1/iV. These 
are informally referred to as cat states in honor of the 
famous Schrodingcr cat. A particular set of cat states 
is generated by passing a coherent state through a 
Kerr medium, i.e. applying the Hamiltonian evolu- 
tion corresponding to 

H K = Ca-N + 7 kN 2 . 

Namely, \ip) = e~ ltKHK is a cat state of type (13) 
with N components if txTK = tt/N . Unlike coherent 
states, cat states are very quickly destroyed by typ- 
ical perturbations. Therefore, the in practice slow 
procedure of generating a cat state through a Kerr 
medium is poorly usable, and more subtle schemes 
must be devised to stabilize them. 

In [18, 17], by tailoring the initial atom state (u) 
and imposing a simple time-varing profile of the inter- 
action parameter (5(t)), we were able to propose an 
engineered reservoir such that one atom-field inter- 
action, lasting a time T, applies the following Kraus 
map to the field (modulo some reasonable approxi- 
mations; see [17] for details): 



P 
M k 

Mi 
Mi 



MxpM\ + M 2 pM\ with 
e - ihN M k e ihN for A; = 1,2: 



(14) 



cos(f) cos(^) +sin(§) 



sin(|) cos(- 



t)-cos(f)a 



N 

sin(^H) 



N 



/(0)N 



Here the photon annihilation operator is defined as 
a = Y^n=n V™\ n ~^)( n \ an d satisfies a^a = N. Pa- 
rameter u € [0, 2ir] reflects initial atomic state choice, 



while 9 and <f) can be varied in [0, 2ir] by properly 
tailoring 6(t) during the interaction. The trick for 
"simulating" the Kerr-Hamiltonian-like factor is 
to play with the non-commutation of evolutions as- 
sociated to the different values that S(t) takes over 
time. 

Lemma V.l: The dynamics of p under the Kraus 
map associated to M\,M2, corresponds to dynamics 
for p = e t ' lN ( oe _l ' lN according to the Kraus map asso- 
ciated to Mi, M 2 - That is, the state trajectory p(t) is 
equivalent to the state trajectory p(t) viewed through 
a Kerr medium. □ 

Thus if Mi , Mi stabilize a state poo close to a 
coherent one - this is a not too exotic task - then the 
actual interaction can stabilize a state poo close to 
a cat state. In the following, for simplicity we focus 
on the case <fi = tt corresponding to a 2-component 
cat; a similar analysis can be carried out for a more 
general case. 

If u, 9 are small, the Kraus maps for p and p are 
close to the identity and can be viewed as discretiza- 
tions of Lindblad differential equations. One readily 
computes: 

f t p = [/?at - /?t a , p] - |(Np + pN - 2apat) , (15) 

with f3 dt = u9/A and k dt = 9 2 /4. To this we can add 
Lindblad terms corresponding to disturbances. The 
dominant disturbance is spontaneous photon loss, 
corresponding to a term — 4f (Np + pN — 2a/oa^) in 
the original frame. Transforming that term to the 
frame of p and adding it to (15) yields the dynamics: 



TtP 



Pl&t - a, p] - § (Np + pN- 2apa t ) 
-f (Np + pN- 2e" rN a / 5a t e-" rN ) . (16) 
/3[a f - a, p] - ^(Np + pN - 2apa t ) 
-K c {apaJ - e" N - 



-apaV" N ) 



The last term induces a structural difference between 
equations (15) and (16); due to the change of vari- 
ables, it implies non-local interference since an oper- 
ator is compared to its rotation by ir in electromag- 
netic phase space. 
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5.3 Convergence analysis 

We analyze the convergence properties of p under 
(16); the corresponding properties for p are identi- 
cal, through the unitary change of frame e~ l ' iN . 

Lemma V.2: The dynamics (16) has an equilibrium 
of the form 



jji{z) \z)(z\ dz 



(17) 



where \ z) is a coherent state with z G R, a c = 2/3/ (n- 
k c ), and 



p(z) = po 



(« 



2 2\ (2k c )/(k + K c ) 

a 2 _ z 2 a c e z 



a c - z 



(18) 



with po a normalization constant ensuring 
f° c ac H{z)dz = l. 

Proof: One easily checks that (16) leaves invari- 
ant the set of states of the form (17) with arbitrary 
p(z). Plugging this form into (16) and looking for a 
stationary solution p(z), a, projection on all real co- 
herent states followed by a resolution by parts of the 
resulting ordinary differential equation (with nonlon- 
cal terms!) yields the result; see Appendix B of [17] 
for details. □ 

In any case, p{— a c ) = 0. At the limit k c — > 0, the 
distribution p{z) converges to a Dirac distribution at 
z = a c , implying: 

Corollary V.3: The dynamics (16) without distur- 
bance, i.e. with k c = 0, has an equilibrium of the form 
Poo = \otc){a c \, for which p = p^ is a two-component 
cat state of the form (13), more precisely: 



(|d c ) + i\-d c )) (\d c ) + i\-d c )) 



with \d c \ = |a c |. 

(This result is more easily obtained from (15) di- 
rectly.) □ 

For small k c , we have lim z _>. Qc p{z) = +oo so 
Poo is close to the coherent state |a c )(a c |. We thus 
approach our goal. However, it remains to assess 
how/whether an arbitrary initial state converges 
towards poo. For this we use the Corollary IV. 3 



developed above. 

For starters, note that p^ has full rank as long 
as k c 7^ 0. Indeed, the set of coherent states {|a) : 
a € [a,b] C K} makes an overcomplcte basis of our 
Hilbert space, for any b > a. An annoying point is the 
infinite-dimensional setting. Since n represents the 
energy of Fock state \n) in units of hui, physical argu- 
ments make it acceptable to limit ourselves to a finite- 
dimensional Hilbert space, spanned by |0), |1), |n): 
nobody physically expects the photon loss to trans- 
form the stable, coherent-state-stabilizing dynamics, 
into one where on average a significant fraction of the 
state is lost to infinitely high energies. Within this 
approximation, we can then apply Corollary IV. 3 di- 
rectly to get the following. 

Proposition V.4: Take n > and de- 

note by P fi the orthogonal projection onto Wn = 
span{\0),\l},...,\n}}. Set N ft = P S NP S and a n = 
P fi aP fi the truncation to o/N and a. The trun- 
cation of (16) to Hn reads: 

= 

dtrn 

,9 [at - an, p n ] - '^C^nPn + p fi N fi - 2a fi p fi at) 



Kc{&-nPn a " 



e An Pn A fi e 



(19) 



where the support of pn is included in T-L n - When (19) 
admits a full rank steady state p n ,oo, the Lyapunov 
function VBures(Pn) applied to (19) is strict, showing 
that the truncated system globally converges towards 

Pn,oo • 

Proof: From Corollary IV. 3, we have that 
VBures(Pn) is a Lyapunov function for any Lindblad 
equation. It remains to examine the set where 

it ^Burcs = -Tr( [G, a fi ] p R)0O [G, a^] f ) 

- Tr( [G, e" Nfi a fi ] p n ,oo [G, e^a*]* ) = . 

With Pn,ao positive definite, this occurs only when 
G commutes with & n and with e l7vNii & n ; in fact one 
of them will be sufficient for our proof. Since G is 
self-adjoint, we have the conditions Gag = a^G and 
also Gat = ajjG, which together yield N fi G = GN n 
since N s = ata fi . Thus G is diagonal in basis 
{|0), |1), \n)}, and as it commutes with a s it must 
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in fact be a multiple of identity, G = AI,-, . Plugging 
this into (11) and taking the trace shows that neces- 
sarily A = 0, thus G — which implies p fl = Pn l00 . 
□ 

Thanks to the tools developed in Section IV, the 
proof above makes notably little use of system- 
specific properties. It can in fact be repeated ver- 
batim to show the following. 

Theorem V.5: Consider any Lindblad-Kossakowski 
system (3) which features one decoherence term Lk 
proportional to any finite dimensional truncation a„ 
of a (defined in Proposition V.4) and has a full rank 
equilibrium p^. Then Vsuresip) defined in Corollary 
IV. 3 is a strict Lyapunov function for this system, so 
the latter globally converges to p^. 

Note that to have the same result with the DFS cri- 
terion [3, 20, 5], it is necessary to assume in addition 
that the subspace generated by all Lk is Hcrmitian. 
This would be obtained e.g. with a thermal bath at 
nonzero temperature. It seems worthwhile that The- 
orem V.5 gives a formal result also for environments 
at zero temperature, which is the ideal situation for 
physical experiments. 

6 Conclusion 

On the basis of a set of contractive metrics for 
Kraus maps, we have systematically built a set 
of Lyapunov functions which decrease under any 
Lindblad-Kossakowski evolution. In particular, 
unlike for many existing and popular convergence 
criteria such as fidelity, we obtain practical expres- 
sions for the rate of decay, allowing an efficient global 
investigation of mixed stationary states of the evolu- 
tion. A major restriction of our approach is that it 
assumes an equilibrium of full rank. At the opposite 
side, pure state equilibria can usually be treated 
reasonably with fidelity-like indicators. An open 
question is thus how to efficiently treat systems with 
steady states of partial rank. We have illustrated our 
tool by proving convergence of an engineered reser- 
voir scheme that stabilizes 'Schrodinger cat' states 
of an electromagnetic field mode. For this we have 
introduced an approximation of finite dimensional 



Hilbert space. Infinite-dimensional Hilbert spaces, 
related to nonlocal PDE formulations, are typical of 
quantum systems like the basic harmonic oscillator. 
Systematic tools to alleviate the problems related 
to infinite dimension would thus be a welcome 
completion of this work. 
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